Microscopic Calculation of Fusion Cross-Sections 
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A microscopic calculation of cross sections for fusion of oxygen isotopes ie O, 22 and 24 is pre- 
sented. Fermionic Molecular Dynamics wave functions are used to describe the fully antisymmetrized 
and angular momentum projected nucleus-nucleus system. The same effective nucleon-nucleon in- 
teraction is used to determine the ground state properties of the nuclei as well as the nucleus-nucleus 
interaction. From the microscopic many-body wave function the corresponding wave function for 
the relative motion of two point-like nuclei is derived by a method proposed by Friedrich which leads 
to a local effective nucleus-nucleus potential. Finally the Schrodinger equation with incoming wave 
boundary conditions is solved to obtain the penetration factors for the different partial waves. With 
these the S-factor for the fusion process is calculated. A good agreement with experimental data is 
obtained for the 16 0- 16 cross section. Much enhanced cross sections are found for the neutron-rich 
oxygen isotopes. 

PACS numbers: 21.60.-n,21.60.Gx,25.60.Pj,26.50.+x 



I. INTRODUCTION 

Heavy-ion fusion reactions play important roles in var- 
ious astrophysical objects. For example, reactions in- 
duced by the fusion of two 16 O nuclides are the main en- 
ergy source of hydrostatic oxygen burning [l|, Q • On the 
other hand pycnonuclear fusion reactions of neutron-rich 
nuclides are expected to occur in the crust of accreting 
neutron stars in binary systems, when the ashes of x-ray 
bursts are pressed to high densities and a series of elec- 
tron captures has transformed the matter neutron-rich 
[3, [1] . Electron captures and pycnonuclear fusion gener- 
ate energies which can affect the dynamics and evolution 
of the neutron star crust Q . 

To determine the relevant fusion cross sections is quite 
challenging. While for the 16 0+ 16 reaction data ex- 
ist at energies below the Coulomb barrier d, 0, H, Q, 
it is yet impossible to determine the cross sections at 
the most effective energies of hydrostatic oxygen burning 
(£k2 MeV) . Hence it has been the standard procedure 
to extrapolate the measured cross sections down into the 
energy region of astrophysical interest usually based on 
phenomenological barrier penetration models (ic| . Fur- 
thermore, a calculation of the subbarrier 16 0+ 16 fusion 
cross sections performed on the basis of the Hartree-Fock 
model agrees with the data quite well . The situation 
is significantly more uncertain for the determination of 
the fusion cross sections for the neutron-rich nuclei. Here 
no data exist and the cross sections have to be derived 
solely from theoretical models (e.g. 12]). Furthermore 
neutron-rich nuclei develop neutron densities with long 
tails which lead to strongly enhanced fusion probabili- 
ties, but are very difficult to systematically describe by 
parametrized potentials [l2j. 

Obviously it is desirable to base the required extrapo- 
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lations on theoretical models which are as reliably as pos- 
sible; i.e. which are accounting for the relevant degrees of 
freedom of the reaction. In recent years powerful many- 
body approaches have been developed (Green's Function 
Monte Carlo [Hj], No-Core Shell Model [Hj], Fermionic 
Molecular Dynamics (FMD) [IE (MODI) which allow now 
a quantitatively correct description of light nuclei based 
on the solution of the many-body problem for realistic 
nucleon-nucleon interactions. The next step should be 
to take advantage of these advances in nuclear structure 
physics for a more realistic and reliable description of low- 
energy reaction problems. First attempts have been pre- 
sented in [let fl3 | where microscopic structure informa- 
tion derived from the GFMC or no-core shell model has 
been incorporated into potential model descriptions of 
low-energy d+d and p+ 7 Be fusion. Here we will present 
a different approach, based on the FMD model, which al- 
lows for a consistent description of bound, resonant and 
scattering states from the same many-body Hamiltonian 
derived from a realistic nucleon-nucleon interaction. In 
particular, our model describes consistently the internal 
structure and properties of the asymptotic reaction part- 
ners and the relevant degrees of freedom of the collision 
process. This is particularly important for astrophysi- 
cally relevant reactions including nuclei far from stabil- 
ity, as their description is usually not constrained by data 
yet, but it is expected to depend sensitively on nuclear 
structure subtelities like the formation of neutron halos 
or skins. While the FMD approach to nuclear reactions 
is quite general and can be applied to scattering, cap- 
ture and transfer reactions, we will here present a first 
application which deals with subbarrier fusion of oxygen 
isotopes. 

We introduce the FMD many-body model in Sec. [Til 
The calculation of the microscopic nucleus- nucleus po- 
tential and the transformation into an effective potential 
for point-like nuclei is performed in the language of the 
well known Generator Coordinate Method (GCM) and 
Resonating Group Method (RGM) and is presented in 
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E b I MeV (Experiment) -127.62 -162.03 -168.48 
E b I MeV (FMD) -120.07 -153.88 -158.75 

rwttor / frn (Experiment) 2.54(2) 2.88(6) 3.19(13) 
r-mattcr / fin (FMD) 2.50 2.80 2.95 

TABLE I: Ground state energies and matter radii calculated 
with the FMD wave functions. The experimental matter radii 
are taken from 12311. 



Sec. [TTT] and [IV] In Sec. [V] the fusion cross sections are 
calculated by solving the Schrodinger equation using in- 
coming wave boundary conditions. In Sec. IVII wc discuss 
adiabatic effects before summarizing in Sec. IVII1 




II. MANY-BODY MODEL 

In the Fermionic Molecular Dynamics (l5l . [l6| model 
the many-body states are described by Slater determi- 
nants 

\ip) = -4||«i>®...® |gyt)j (l) 

using Gaussian wave packets for the single-particle states 

(x\q) ^pI- ^^ } ® |x T ,X X ) ® |0 ■ (2) 

The complex parameter b encodes the mean position and 
the mean momentum of the wave-packet. In this paper 
the width parameter a is chosen to be identical for all 
single-particle states. The spins can orient freely, the 
isospin is fixed to describe either a proton or a neutron. 
The ground state wave function is obtained by minimiz- 
ing the intrinsic energy with respect to the single-particle 
parameters a, bi and Xi- For 16 0, 22 and 24 we obtain 
configurations that are identical to closed-shell harmonic 
oscillator configurations with a fully occupied 0d 5 / 2 shell 
in 22 and a fully occupied lsi/2 shell in 24 0. 

We use a nucleon-nucleon interaction that is derived 
from the Argonne V18 interaction by explicitly introduc- 
ing the short-range central and tensor correlations into 
the many-body system by means of the Unitary Correla- 
tion Operator method [20 , l2l| in a two-body approxima- 
tion. The correlated interaction has been used success- 
fully in ab initio calculations using many-body pertur- 
bation theory and the no-core shell model [22|. As the 
simple FMD wave functions are not able to describe the 
medium- to long-range tensor correlations we use a phe- 
nomenological modification of the two-body interaction 
as introduced in Here an additional momentum- 

dependent two-body part is used and fixed to reproduce 
the binding energies and radii of 4 He, 16 and 40 Ca. 
An additional spin-orbit term is tuned to the binding- 
energies of 24 0, 34 Si and 48 Ca. 



FIG. 1: Matter densities for ie O, 22 and 24 0. Finite size of 
protons and neutrons is taken into account. 
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FIG. 2: Charge densities for 16 0, 22 and 24 0. Finite size of 
protons and neutrons is taken into account. 



The calculated binding energies and matter radii are 
shown in Tab.[U The ground state densities are shown in 
Figs. [U [U An improved description of the ground states 
is possible in a variation after projection approach. But 
as the intrinsic states would then no longer be angular 
momentum eigenstates, the calculation of the nucleus- 
nucleus potential would no longer be feasible. As the 
density distribution remains almost unchanged we expect 
no significant effect on the nucleus-nucleus potentials and 
therefore on the fusion cross sections. 
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III. WAVE FUNCTIONS FOR RELATIVE 
MOTION 

A. GCM wave functions 

To calculate the microscopic nucleus-nucleus energy 
surfaces we make use of the fact that the FMD wave 
functions can be easily moved and rotated. We construct 
Generator Coordinate Method (GCM) wave functions us- 
ing the normalized FMD ground state wave functions 
and the distance R = Re z between the origin of the slater 
determinants as the generator coordinate. R is assuming 
values from 0.5 fm up to 15 fm in steps of 0.5 fm 



\V(R)) = A{\M-R/2)) ® |^(+i?/2))} 



(3) 



We have to project these states on angular momentum 
with the projection operator [241 ] 
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(4) 



il = (a,[3,j) denoting the three Euler angles. As the 
intrinsic spin of the nuclei is + the total angular mo- 
mentum J is given by the orbital angular momentum 
L of the relative motion of the two nuclei. As the two 
nuclei are identical we will only encounter even angular 
momenta and positive parity states. Because of the ax- 
ial symmetry of the system only the integration over the 
azimuthal angle j3 has to be carried out explicitly which 
we do using a Gauss-Legendre scheme with 50 points. 
Because of the property 
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we can calculate the projected GCM Hamiltonian and 
overlap functions as 



H L {R, R') 
and 



——^(i^) |(ff - T cm )P^{R'e z )) 

(6) 



Air 

N L (R, R') = ^—^(^(Re z ) \P^(R'e z )) . (7) 



The calculated GCM energy surfaces 



E L {R) 



H L (R,R) 
N L (R,R) 



(Ebi + Eb2) 



(8) 



where Ebi and Eb2 are the ground state energies of the 
fragment nuclei, are shown for L = in Fig. [31 The GCM 
energy does not approach the Coulomb potential at large 
separations as the kinetic energy due to localization of 
the nuclei is still included. At small separations we can 
observe the effects of antisymmetrization and saturation 
that are missing in a double-folding potential. 




FIG. 3: GCM energy surfaces for L — 0. The energy surfaces 
are different even asymptotically because of different localiza- 
tion energies in the relative motion. 



We can fully describe the system of the two nuclei using 
the GCM basis states Eq. ([3j) 



*lm) = J dRR 2 P MQ \t>(Rel))f(R) 
*I>j&o|*(fli£)>/<- 



(9) 



But the GCM "wave function" f(R) can not be inter- 
preted as a wave function for the relative motion of two 
nuclei. This is because of the non-orthogonality of the 
basis states. Friedrich [25| proposed a two-step approach 
for dealing with the non-orthogonality. In a first step 
a transformation from the GCM into the Resonating 
Group Method (RGM) basis is performed. This elim- 
inates the non-orthogonality that is caused by the use 
of Slater determinants in the GCM. In a second step 
the non-orthogonality due to the antisymmetrization be- 
tween the two nuclei at short distances is removed. 



B. RGM wave functions 

In the RGM the wave function for the relative motion 
of the two nuclei is given using basis states ^(r 1 )) 

(P,ti,b\*(i f )) = A{6(f f -fl(l>(S l )4>fa)} . (10) 

Here (f>(£) are the intrinsic wave functions of the nuclei. 
If all the single-particle states in the GCM Slater 
determinants have the same width parameter a the GCM 
basis states can be expressed with the RGM basis states 
as MM 



d 3 rr(i?-r)|$(r))|$ cm ) (11) 
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Collective wave functions 



/ a \ 3 / 4 
r(fl-f)=( — ) exp 
\7ra/ 



2 a 



M = 



AiA 2 
A x + A 2 
(12) 



and the center-of-mass wave function 



<X|$ cm ) 



3/4 



exp —A 



2a 



A = A x + A 2 



(13) 

We introduce the angular momentum projected basis 
states 



The RGM norm kernel 

n L {r,r') = ($LAi(r) \3>LAl(r')) 
asymptotically behaves as 

/, r,r'^oo 1 + <5l2 ( — 1) 



n L (r,r ) 



<5(r — r ) 



(14) 



(15) 



(16) 



Asymptotically the RGM wave function (|2lj) possesses 
the properties of a "proper" wave function and can for 
example be used for implementing boundary conditions. 
At short distances the RGM norm kernel Eq. (fT5|) de- 
viates from the asymptotic behavior because of the an- 
tisymmetrization between the two nuclei. To obtain a 
wave function that can be identified as a wave function 
for the relative motion of two point-like nuclei we have 
to transform the RGM basis states |$lm(?")) new 
basis states |$lm(^)) that fulfill 

($LM(r)\$ LM (r')) = 1 + Sl *t 1)I < r - r ') ■ ( 22 ) 
This is achieved by 

\*LM(r)) = J dr>r' 2 \$ LM (r>)) nl 1/2 (r',r) . (23) 

For the evaluation of the norm kernel we use its spectral 
representation 



n L (r,r') = J2x {a) (r) ^ X (a \r') 



(24) 



In the case of identical nuclei only even L are allowed but 
an additional factor of two appears in the normalization 
(exchanging all the nucleons between the two nuclei does 
not change the wave function). 

The angular momentum projected GCM basis state 
can be expressed by the angular momentum projected 
RGM basis state as 

Pk I0 \V(Re z )) = J drr 2 u L {R;r)\<5> LM {r))\§ cra ) (17) 



with 



(18) 



and 



/ u \ 3 / 4 
T L (R;r) =47r(-£-) exp 
\Tra/ 



3/4 / R 2 + r 2 \ ( Rr 
2a I \ a 



19) 



where il are the modified spherical Bessel functions [2£ 

The system of two nuclei (JS]) can now be described 
using RGM basis states 

|*lm)= /drrV GM M|$LAf(r))|$ cm ) (20) 



(21) 



with the RGM wave function 



RGM 



(r) = ^u L {R t ;r)f l 



where (j,( a > are the eigenvalues and x are the (real) 
eigenvectors of the norm kernel. We perform the diago- 
nalization numerically using the GCM basis states. For 
that we expand the eigenvectors 



X (o) (r)=5^u £ (JJ i ;r) 



,(«) 



and transform the integral integration 

dr'r' 2 n L (r,r')x (a) (r') = fJ- (a) X {a) (r) 
into the generalized eigenvalue problem 



j 



(25) 



(26) 



(27) 



with 



drr / dr' r' UL(Ri;r)n L (r,r')u L (Rj;r') 



(^{RA^P^iR^)) 



and 
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dr r 2 u L {Rt\ r)u L (Rj;r) . 



(28) 



(29) 



We can now express the GCM basis state as 



P& \*(R l e' a )) = I drr 2 u L (Ri;r)\$ LM (r))\$> cm ) (30) 
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FIG. 4: The wave functions rv,L(R;r) (dashed lines) and 
ru L {R;r) (solid lines) for 16 0- 16 0, L = at R = 4.5 fm 
and R — 11 fm. ul{R,t) is suppressed at small distances due 
to antisymmetrization. 



with 



u L (i?,;r) = / dr' r' 2 n^ 2 (r,r') u L {Ri;r') 



k I \ a / 



Vl ■ 

(31) 



The effects of antisymmetrization are illustrated in 
Fig. [4] where ui,{R;r) and ui,{R)r) are shown for small 
and large separations. At large separations ul and co- 
incide while at R = 4.5 fm, where there is a large overlap 
between the surfaces of the nuclei, u is substantially re- 
duced and shifted outwards due to antisymmetrization. 
Completely Pauli forbidden states correspond to eigen- 
values n = 0. For large separations the antisymmetriza- 
tion between the clusters has no effect and the eigenvalues 
fi approach 1. 

For a state described in the GCM world as 



\* LM )=Y, P M^{Rie z ))fi 



(32) 



the corresponding proper wave function for the relative 
motion is given by 



<Pi<(r) = ^2u L (Ri;r)fi 



(33) 



Friedrich [25j and fit a local potential to the GCM matrix 
elements. 

When comparing the GCM overlap matrix elements 
with the corresponding expression in the collective world 

(34) 

we observe numerical deviations for Ri < 3.5 fm where 
Ul(R4',t) is almost completely suppressed. As we will 
solve the Schrodinger equation with incoming wave 
boundary conditions, we will need the solution of the 
Schrodinger equation only from the minimum in the po- 
tential surface outwards to calculate the barrier penetra- 
tion. We therefore fit the Hamiltonian matrix elements 
only for Ri > 4.0 fm. Our ansatz for the effective Hamil- 
tonian in the collective world is given by 



n 2 

Picl 



L(L+1) 

2fj,M N ^ 2[iM N r 2 



+ V e L s (r) + V c (r) + E bl +E b2 

where the Coulomb interaction is that of two homoge- 
nously charged spheres (r c is taken to be the charge ra- 
dius) 

\z x Z 2 e 2 \ ;r>2r c 
^ZiZ 2 e z (3 - 77)27- \r<2r c 

The effective potential Vjfair) is parameterized as a 
sum of ten Gaussians. We minimize the deviation from 
the GCM matrix elements to fulfill the relation 

{^(R t e z ) \(H ~ T cm )P^(R t e- z )) = 

drr 2 u L (Ri;r)H^{r)u L (Ri;r) (37) 



as well as possible. We determine Vjg(r) independently 
for L = 0, 2, 4 and obtain almost identical fits. Therefore 
a L-independent V c g(r) is chosen and fitted simultane- 
ously to the L — 0, 2, 4 matrix elements. The deviations 
between the collective and GCM matrix elements are less 
than 50 keV. 

The fitted L = effective potentials including the 
Coulomb interaction are shown in Fig. [5] As could al- 
ready be deduced from the GCM energy surfaces the 
neutron rich isotopes feature lowered barriers at larger 
distances due to the attractive nuclear interaction in the 
tails of the neutron distribution. 



V. FUSION CROSS SECTIONS 



IV. NUCLEUS-NUCLEUS POTENTIALS 

In principle the transformation procedure defines an 
energy independent effective interaction but it has a com- 
plicated non-local structure. We will therefore follow 



With the effective potential derived in the last section 
the two-body Schrodinger equation 



L(L + 1) 
2(iM N ' 2fiM N r 2 



Prel 



+ V e s(r) + Vc(r) 



Ecp L (r) 
(38) 
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FIG. 5: Effective Nucleus-Nucleus Potentials [L = 0) derived 
from the FMD energy surfaces. The barrier gets lower and 
shallower with increased neutron number. 



is solved with Incoming Wave Boundary Conditions 
(IWBC) [III. No imaginary part of the potential is 
needed in this approach. We assume that the nuclei will 
fuse when the minimum in the potential surface has been 
reached. The IWBC can be formulated as 



?fci(r m i n ) 



where k^ir) is the local wave number 



(39) 



k L (r) 



' 2fj,M N E - 



L{L + 1) 
2fiM N r 2 



Kff(r) - V c (r] 



The Schrodinger equation is integrated numerically 
starting from r m j n outwards and matched to the Coulomb 
scattering solution 



C-(F L (kr)-iG L (kr)) + 
C+(F L (kr)+iG L (kr)) 



with 



k = yj2[iM N E . 



(41) 



(42) 



The matching point r m has to be outside the range of 
the nuclear potential V e g(r). The penetration factor Pl 
is given by the ratio of the incoming flux at r m j n to the 
incoming Coulomb flux 

Pl 



kL(r min) 



k\C- 



(43) 



For very low energies (E < 4 MeV) the magnitude of 
the transmission coefficient gets very small and the nu- 
merical matching is no longer possible. For these energies 
we use the approximation [29j 



k L{r min) |</?L(?"min)| 2 



AFl{kr v 



k L (r m ) \fL{r m )\ 2 Fl(kr m ) + G\{kr m ) 



(44) 



With the penetration factors we immediately obtain 
the fusion cross section 



< E ) = p E( J + S ^(-l) L )(2L + l)P L {E) . (45) 

L=0 

As we are dealing with the fusion of identical nuclei 
only even angular momenta L contribute to the cross 
section. 

Instead of the cross section we plot in Fig. [6] the astro- 
physical S-factor where the tunneling through a point- 
charge s-wave Coulomb barrier has been taken out of the 
cross section 



S{E) = a(E)Ee 



2-KT) 



(46) 



For 16 O the agreement with the available data is satis- 
factory, taking the fact into account that no parameters 
have been adjusted. The observed edge in the S-factor 
is related to the barrier height (see Fig. [5]) . In the case 
of 24 O no fusion occurs below E cm = 2 MeV because the 
potential calculated with frozen configurations does not 
drop below zero inside the Coulomb barrier. In the fol- 
lowing section we will discuss the adiabatic energy surface 
for which the S-factor will smoothly continue to E cm = 0. 

Comparing the neutron rich oxygen isotopes with 1 O 
we find fusion cross sections that are enhanced by 6-8 
orders of magnitude. This is caused by the reduction 
of the height and width of the Coulomb barriers for the 
neutron rich 22 O and 24 0. 

In the last years a hindrance of the fusion process at 
sub-barrier energies compared to coupled channel calcu- 
lations has been found for medium- heavy nuclei ( (30j and 
references therein). In [3l[ this effect was related to a 
shallower potential due to the saturation properties of 
nuclear matter that become important for large overlap 
of the reacting nuclei. In a recent paper [32[ a drop-off 
of the S-factor in the 16 0- 16 fusion reaction has been 
predicted for energies below 7 MeV due to this fusion 
hindrance effect. Although our wave functions are fully 
antisymmetrized products of frozen ground states and 
therefore saturation effects are included explicitly we do 
not find such a decrease in the S-factor. The nucleus- 
nucleus potential while featuring only a shallow minimum 
still allows for an increase in the S-factor for energies well 
below the barrier. 



VI. 



ADIABATIC EFFECTS 



Up to now we have discussed the fusion process in 
a single- channel approximation where the nuclei stay in 
their ground states. To describe effects beyond the single- 
channel approximation we could in principle extend our 
discussion to a coupled-channel approach by including 
excited states and trying to derive the coupling poten- 
tials. 
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FIG. 6: S-factor for the fusion of oxygen isotopes. Experi- 
mental data I1,0,HH for ls O are included. The S-factor is 
calculated with the effective potential derived for the frozen 
configurations (thin lines) and with the modified potentials 
for the adiabatic configurations (thick lines). 
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FIG. 7: GCM energy surfaces for the 24 0- 24 system 
calculated with frozen configurations (solid line) and with 
FMD configurations minimized under quadrupole constraints 
(dots). The intrinsic quadrupole deformation is mapped on 
the separation R between the nuclei. The effective nucleus- 
nucleus potential derived from the frozen GCM configs is 
given by the dashed line. 

Another point of view is to switch from the diabatic 
picture with frozen configurations to an adiabatic ap- 
proach where we allow for distortions of the wave func- 
tions. We generate such wave functions by minimizing 
the expectation value of the energy for the 2A-nucleon 
system with respect to the parameters of all the FMD 



single particle states under constraints on the quadrupole 
deformation of the total intrinsic wave function. For 
these configurations we define the separation distance 
to be that of the frozen configuration with the same 
quadrupole deformation. The adiabatic GCM energy 
surface calculated for 24 can be seen in Fig. [7] We 
see almost no effect in the region of the Coulomb barrier 
but observe less repulsion at smaller separations. The 
corresponding adiabatic wave functions can no longer be 
expressed with RGM wave functions and therefore the 
derivation of the effective nucleus- nucleus potential as 
done for the frozen configurations is not possible any 
more. To estimate the effects of the adiabatic configura- 
tions we modify the nucleus- nucleus potential as shown 
in Fig. [7] to reflect the changes seen in the GCM energy 
surface. These modifications change the low energy be- 
havior especially in the 24 0- 24 case where the nucleus- 
nucleus potential for the frozen configurations has a min- 
imum above zero energy. The modified potential goes 
below zero and the S-factor now smoothly continues to 
E cm = 0. A similar but less pronounced effect can be 
seen for 16 0- 16 0, whereas the 22 0- 22 S-factor stays al- 
most unchanged. Variations of the depth of the modified 
potential from — 0.5MeV to — lOMeV change the S-factor 
by less then 50%. 



VII. SUMMARY AND CONCLUSION 

In this paper we studied the fusion cross sections for 
the oxygen isotopes 16 0, 22 and 24 in the Fermionic 
Molecular Dynamics model. The same effective nucleon- 
nucleon interaction is used for the ground states and 
the two-nucleus system. The total wave function is 
antisymmetrized and projected on angular momentum 
eigenstates. In a two-step process an effective nucleus- 
nucleus potential is derived from the GCM energy sur- 
face. The Schrodinger equation is solved with incoming 
wave boundary conditions to obtain the fusion cross sec- 
tions. Effects beyond the frozen ground state approx- 
imation are studied with adiabatic configurations ob- 
tained as FMD many-body states under constraints on 
the quadrupole deformation. Experimental data on 16 O- 
16 are reproduced; a much enhanced cross section for 
the fusion of the neutron-rich oxygen isotopes is found. 
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